Information theory of open fragmenting systems 
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An information theory description of finite systems explicitly evolving in time is presented. We 
impose a MaxEnt variational principle on the Shannon entropy at a given time while the constraints 
are set at a former time. The resulting density matrix contains explicit time odd components in 
the form of collective flows. As a specific application we consider the dynamics of the expansion in 
connection with heavy ion experiments. Lattice gas and classical molecular dynamics simulations 
are shown. 
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I. INTRODUCTION 

The microscopic foundations of thermodynamics are 
well established using the Gibbs hypothesis of statistical 
ensembles maximizing the Shannon entropyQ. At the 
thermodynamic limit, the various Gibbs ensembles con- 
verge to an unique thermodynamic equilibrium. How- 
ever, most of the systems studied in physics do not cor- 
respond to this limit 0. In finite systems the various 
Gibbs ensembles are not equivalent |3| and the physical 
meaning and relevance of these different equilibria has to 
be investigated. 

A common interpretation of a statistical ensemble for a 
finite system is given by the Boltzmann ergodic assump- 
tion. In this interpretation the statistical ensemble repre- 
sents the collection of successive snapshots of a physical 
system evolving in time, and the state variables are iden- 
tified with the conserved observables. This interpretation 
suffers from important drawbacks. First, even for a truly 
ergodic Hamiltonian, a finite time experiment may very 
well achieve ergodicity only on a subspace of the total 
accessible phase spaceQ. Moreover, ergodicity applies 
to confined systems and thus it requires the definition 
of boundary conditions . Then the statistical ensemble 
directly depends on the boundary conditions and we will 
discuss that an exact knowledge of the boundary corre- 
sponds to an infinite information and is therefore hardly 
compatible with the very principles of statistical mechan- 
ics . Finally, the systems experimentally accessible are 
often not confined but freely evolve in the vacuum, as it 
is notably the case for heavy ion collisions. The concept 
of a stationary equilibrium defined by the variables con- 
served by the dynamics in a hypothetical constraining 
box, is not useful for these systems. 

However statistical approaches, expressing the reduc- 
tion of the available information to a limited number of 



collective observables, arc still pertinent to complex sys- 
tems even if the dynamics does not allows at any time 
a total and even exploration of the energy shellpl. In 
this interpretation the MaxEnt postulate has to be inter- 
preted as a minimum information postulate which finds 
its justification in the complexity of the dynamics inde- 
pendent of any time scale 0,0. 

This information theory approach is a very powerful 
extension of the classical Gibbs equilibrium: any arbi- 
trary observable can act as a state variable, and all sta- 
tistical quantities can be unambiguously defined for any 
number of particles |6j. The price to be paid for such a 
generalization is that the constraining variables as well 
as the density matrix continuously evolve in time. The 
time dependence of the process naturally leads to the ap- 
pearance of new time-odd constraints or collective flows. 
In the case of an ideal gas of particles or clusters we 
will show that the ensemble of constraints forms a closed 
algebra and that the information at the initial time is 
sufficient to calculate the exact density matrix at any 
successive time. 



II. STATISTICAL EQUILIBRIA 

When the system is characterized by L observables 
known in average < Ag >= TvDA?, statistical equilib- 
rium corresponds to the maximization of the constrained 
entropy 

S c = -THMog-D -^2\ e <A e > 

e 

where D = J2(n) \ ^ {n) ) P {n) (* (n) | is the density ma- 
trix, and A = {A^} are Lagrange multipliers. Gibbs equi- 
librium is then given by 
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D x = — exp -X.A, 
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where is the associated partition sum. It should be 
noticed that microcanonical thermodynamics can also 
be obtained from the variation of the Shannon entropy in 
the special case of a fixed energy subspace. In this case 
the maximum of the Shannon entropy can be identified 
with the Boltzmann entropy max (S) — log W (E) , where 
W is the total state density with the energy E. In the 
following we shall confine ourselves to the Gibbs formu- 
lation QJ , which is more general than the microcanonical 
ansatz. Indeed the microcanonical density matrix corre- 
sponds to an even occupation of the whole energy shell 
while non ergodic components can be already included 
within the Gibbs formalism through the introduction of 
extra constraints. 



A. Boundary condition problem in finite systems 

The statistical physics formalism recalled above is valid 
for any system size . However, as soon as one Ag con- 
tains differential operators such as a kinetic energy, eq. 
(|T|) is not defined, unless boundary conditions are speci- 
fied. Only at the thermodynamic limit boundary condi- 
tions are irrelevant, as only in this limit surface effects 
are negligible. The definition of any density with a fi- 
nite number of particles requires the definition of a finite 
volume. To this aim, a fictitious container is generally 
introduced The volume and shape of this unphysi- 
cal box has no influence on the thermodynamics of self- 
bound systems, but in the presence of continuum states 
the situation is different. Let us consider the standard 
case of the annulation of the wavefunction on the surface 
S of a containing box V. Introducing the projector, Ps 
, over S and its exterior, the boundary conditions reads 
P s =0 for all microstates (n). Using Pj = P s , we 

can see that this condition imposes an extra constraint 
to the statistical ensemble < Ps >= TrDP s = 0. The 
density matrix then reads 

D Xs = J-exp-ll-6P s (2) 
z \s 

which shows that the thermodynamics of the system 
depends on the whole surface S. For the very same 
global features such as the same average particle density 
or energy, we will have as many different thermodynam- 
ics as boundary conditions. More important, to specify 
the density matrix, the projector P$ has to be exactly 
known and this is in fact impossible. The nature of Ps 
is intrinsically different from the usual global observables 
Ag. Not only it is a many-body operator, but Ps requires 
the exact knowledge of each point of the boundary sur- 
face while no or few parameters are sufficient to define 
the Ag. This infinity of points corresponds to an infinite 
amount of information to be known to define the density 
matrix J3J) • This requirement is in contradiction with the 
statistical mechanics principle of minimum information. 
Thus eq.® is unphysical, and the same is true for the 



standard (N, T, V) or (N, E, V) ensembles when dealing 
with finite unbound unconfined systems. 

B. Incomplete knowledge on the boundaries 

One way to get around the difficulties encountered 
to take into account our incomplete knowledge on the 
boundaries, is to introduce a hierarchy of observables de- 
scribing the size and shape of the matter distribution. 

For example, if only the average system size < R 2 > 
is known, the minimum information principle implies 

bp,p = -J— exp -0 (h + PRoR 2 ) , (3) 

which is akin an isobar canonical ensemble, since the 
additional Lagrange multiplier A^2 imposing the size 
information has the dimension of a pressure when di- 
vided by a typical scale R and by the temperature, 
X R 2 = l3PR . 

A typical application of this concept is the so called 
freeze-out hypothesis in nuclear collisions : at a given 
time to the main evolution (i.e. the main entropy cre- 
ation) is assumed to stop and partitions are supposed 
to be essentially frozen. Typically thermal and chemical 
equilibrium is assumed, meaning that the information at 
to on the energetics and particle numbers is limited to the 
observables < H > and < Nf > for the different species 
/ 013 Freeze-out occurs when the system has expanded 
to a finite size. Then at least one measure of the system's 
compactness should be included . The limited knowledge 
of the system extension leads to a minimum biased den- 
sity matrix given by eq.© |lo| - 

III. MULTIPLE TIME STATISTICAL 
ENSEMBLES 

As soon as one of the constraining observables Ai is not 
a constant of the motion, the statistical ensemble Q is 
not stationary. A single time description may still looks 
appropriate in the freeze out configuration discussed in 
the last section. Indeed in many physical cases one can 
clearly identify a specific time at which the information 
concentrated in a given observable is frozen (i.e. the ob- 
servable expectation value ceases to evolve). However 
this freeze out time is in general fluctuating and differ- 
ent for different observables. For example for the ultra- 
relativistic heavy ion reactions two freeze-out times are 
discussed [9j, one for the chemistry and one for the ther- 
mal agitation. We need therefore to define a statistical 
ensemble constrained by informations coming from dif- 
ferent times. 

Let us now suppose that the different informations on 
the system, < Ae >, are known at different times tf. 
< A( >t,— Trl) (t() Ag. A generalization of the Gibbs 
idea is that at a time t the least biased state of the system 
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is the maximum of the Shannon entropy, considering all 
informations as constraints. 

The maximization of the entropy at time t with the 
various constraints < Ag > tt known at former times tg 
corresponds to the free maximization of 



S c = -Tr If) (t) log D(t)+J2 X deD (tg)j , (4) 

where the Xg are the Lagrange parameters associated 
with all the constraints. This maximization will lead to a 
density matrix which can be considered as a generaliza- 
tion to time dependent processes of the Gibbs ensembles 



Let us consider the case of a deterministic evolution 



d t D = {H,D}, 



(5) 



where { . , . } are Poisson bracket in classical physics and 
commutators divided by ih in quantum physics. The 
minimum biased density matrix is given bv|ll| 



IV. THE DYNAMICS OF THE EXPANSION 

Let us now apply the above formalism to transient un- 
confined systems. Let us concentrate on a scenario often 
encountered experimentally: a finite system of loosely 
interacting particles with a finite extension in an open 
space. We shall assume that at a given freeze out time 
to the system can be modelized as a non interacting 
ensemble of n = 1, ...,iV particles or fragments, and 

a definite value for the mean square radius < R 2 > 

(with R 2 = J2 n ffi) characterizes the ensemble of states. 
Then we have to introduce the constraining observable 

A 2 = R 2 associated with a Lagrange multiplier A . If 
time is not taken into account, the maximum entropy 
solution is given by 
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exp- X.A'[D x (t)], 



(6) 



where the A' £ represent the time evolution of the con- 
straining observables Ag in the Heisenberg representation 
A' e [D (t)] = A' e (At e ) = e- iAt ^A t e iAt ^: 



A' = At 



E 

P =i 



p 



B {p) = {H,B {p - 1] } ■ B {0) =A. 

Eq. @ can be interpreted as the introduction of addi- 

tional constraints B\ and additional Lagrange parame- 

(p) 

ters v) associated with the time evolution of the system 
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-e X p-XrAj2^ p) -B^[D x: ], (7) 



Eq. J7J) is an exact solution of the complete many body 
evolution problem eq.© with a minimum information 
hypothesis on the final time t having made few observa- 
tions < Ag > at previous times tg, which shows the wide 
domain of applicability of information theory. A gener- 
alization of this theory to non deterministic evolutions 
can be found in ref.tLJ. We can see from eq.JZJ) that in 
general an infinite amount of information, i.e. an infinite 
number of Lagrange multipliers are needed if we want to 
follow the system evolution for a long time. However, 
different interesting physical situations exist, for which 
the series can be analytically summed up. In this case, 
a limited information (the knowledge of a small number 
of average observables) will be sufficient to describe the 
whole density matrix at any time, under the unique hy- 
pothesis that the information was finite at a given time. 



Eq. JHJ) is akin to a system of non-interacting particles 
trapped in a harmonic oscillator potential with a string 
constant k = 2Xq/(3. From the partition sum, the EOS 
are easily derived: 



< ft >- 



3m 

T 



< rt >- 



2X ' 



Since Ao-ff 2 is not an external confining potential but 
only a finite size constraint, the minimum biased distri- 
bution © is not stationary. To take into account the 
time evolution, we must introduce additional constrain- 
ing observables 



yz— (fin- 

* — ' m V 



Bp' = - y — ( p n ■ r n + r„ ■ p n 



6(2) _ V- 
R ^m 2 ' 



Since {H,B { ^} = 0, all the other B (p) with p > 2 are 
zero. The density matrix is given by 

A3,Ao (*) = zh: exp En -Peff (*) " A ^ 



2 {^Pn ' ~t~ Tn ' PnJ 7 



(9) 



with 



ef f(t) = P + 2X (t-t o y /m ; u (t) = 2X (t-t )/m 

(10) 

The density matrix 10 can be interpreted as a radi- 
ally expanding ideal gas. Indeed the distribution can be 
written as 
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Dp.x (t) = — — exp y^-fieff (t) 



(p„ - ma (t) f n 



ma 

0.2 0.4 0.6 0.8 1 



A/A, 

0.2 0.4 0.6 0.8 1 



2m 



-Ae// (*) 



(11) 



where a = vq (t) / (3 e f / (t) represents a Hubblian factor 
and the confining Lagrange multiplier is transformed into 



Ae// (t) = A - 



Xo/3m 



2 Peff (t) (3m + 2A (t - i ) S 



(12) 



The term ma (t) r n correcting the momentum in 
cq. (|llf) can be interpreted as a momentum produced by 
a radial velocity a (t) r n . This proportionality of the ve- 
locity with f n shows that the motion is self similar. As 
a consequence, when this collective motion is subtracted 
from the particle momentum, the density matrix 
corresponds at any time to a standard equilibrium © in 
the local rest frame. 

In this case the infinite information which is a priori 
needed to follow the time evolution of the density matrix 

according to eq.|J7|), reduces to the three observables r, 

~p, r ■ p + p ■ r. Indeed these operators form a closed 
Lie algebra, and the exact evolution of (|llf> preserves it 
algebraic structure. 

The description of the time evolution when describ- 
ing unconfincd finite systems has introduced a new phe- 
nomenon: the expansion. One should then consider a 
more general equilibrium of a finite-size expanding finite- 
systems with f3' , a and X' Q as free parameters. Then, if the 
observed minim biased distribution at time t is coming 
from a confined system at time to, the three parameters 
/?', a and X' should be linked to the time to the initial 
tem pera ture /3 _1 and the initial Ao by equations I|1U|) and 

The important consequence is that radial flow is a 
necessary ingredient of any statistical description of un- 
confined finite systems: the static (canonical or micro- 
canonical) Gibbs ansatz in a confining box which is of- 
ten employed misses this crucial point. On the other 
hand, if a radial flow is observed in the experimental data, 
the formalism we have developed allows to associate this 
flow observation to a distribution at a former time when 
flow was absent. This initial distribution corresponds to 
a standard static Gibbs equilibrium in a confining har- 
monic potential, i.e. to an isobar ensemble. 



V. NUMERICAL SIMULATIONS 

As we have already mentioned in section IIVI in the 
hypothesis of negligible interaction between the system's 
constituents the expansion is self-similar, implying that 
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FIG. 1: Upper left: effective pressure (full line) and 
average volume normalized to the ground state vol- 
ume Vo = A (dashed line) as a function of the col- 
lective radial flow. Upper right and lower part: frag- 
ment size distributions in the expanding Lattice Gas 
model. Distributions without flow (full black lines) 
are compared with distributions with 11% (dashed 
black), 43% (full grey) and 150% (dashed grey) con- 
tribution of radial flow at the same temperature (up- 
per right) and at the same total energy (lower part). 



the situation is equivalent to a standard Gibbs equilib- 
rium in the local rest frame. In the expanding ensemble 
the total average kinetic energy per particle is simply the 
sum of the thermal energy (eth) = 3/(2/3) and the radial 
flow (e/i) = ma 2 {r 2 )/2. 

This scenario is often invoked in the literature Q 
to justify the treatment of flow as a collective ra- 
dial velocity superimposed on thermal motion; however 
eq. (|llfl contains also an additional term oc r 2 which cor- 
responds to an outgoing pressure . The phase diagram 
and fragment observables are therefore expected to be 
modified by the presence of flow even in the self simi- 
lar approximation. To quantify this statement, we have 
performed calculations in the Lattice Gas Model^2]> and 
the results are shown in Figure ^ 

The effective pressure A e // as well as the associated 
average volume are shown in the upper left part of fig- 
ure n as a function of the collective radial flow for a 
given pressure and temperature . The Lagrange param- 
eter A e // being a decreasing function of a, the critical 
point is moved towards higher pressures in the pres- 
ence of flow ^tJ- However one can see that the ef- 
fect is very small up to ma ~ .6 (which corresponds 
to Ap 2 /p 2 = (efi)/{e t h) ~ 67% flow contribution). In 
this regime the cluster size distributions displayed in the 
upper right part of figure ^ are on ly slightly affected. 
On the other side if collective flow overcomes a threshold 
value Ap 2 /p 2 « 100% the average volume shows an ex- 
ponential increase and the outgoing flow pressure leads 
to a complete fragmentation of the system (dashed grey 
line in the lower part of fignj. We can also observe that 
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an oriented motion is systematically less effective than a 
random one to break up the system. This is shown in the 
lower part of figure Q which compares for a given A dis- 
tributions with and without radial flow at the same total 
deposited energy: for any value of radial flow equilibrium 
in the standard microcanonical ensemble corresponds to 
more fragmented configurations. 
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time (t ) time (t ) 

FIG. 2: Time evolution of < R 2 >, < P 2 > and 
< R- P > for an initially constrained Lennard Jones sys- 
tem of 147 particles freely expanding in the vacuum, at 
different total energies. Lower right: expansion dynam- 
ics (symbols) compared to the prediction of eq.@ (lines). 

The above results are relevant for the experimental sit- 
uation if and only if the inter-fragment interactions can 
be neglected when the system is still compact enough to 
bear pertinent information on the phase diagram. In- 
deed only in this case the series (J7J can be analytically 
summed up and the expansion dynamics can be reduced 
to a self similar flow 11]. The validity of the ideal gas ap- 
proximation eq.@ for the expansion dynamics is tested 



in Figure E'14] in the framework of classical molecular 
dynamics [15|. A Lennard Jones system is initially con- 
fined in a small volume and successively freely expanding 
in the vacuum. We can see that after a first phase of the 
order of « 10 Lennard Jones time units, where interpar- 
ticle interactions cannot be neglected, the time evolution 
predicted by eq. © is remarkably fulfilled for all total en- 
ergies. This result is due to the fact that the system's size 
and dynamics are dominated by the free particles, while 
deviations from a self similar flow can be seen if the analy- 
sis is restricted to bound particles^. We expect eq.® 
to describe the system evolution even better if the de- 
grees of freedom n are changed from particles to clusters, 
as suggested by the Fisher model of condensation . 



VI. CONCLUSIONS 

In this paper we have introduced an extension of Gibbs 
ensembles to time dependent constraints . This formal- 
ism gives a statistical description of a system observed 
at a time at which the entropy has not reached its sat- 
urating value yet, as it may be the case in intermediate 
energy heavy ion reactions |l3[ . Another physical applica- 
tion concerns systems for which the relevant observables 
pertain to different times, as in high energy nuclear col- 
lisions 0. 

Our most important result is that any statistical de- 
scription of a finite unbound system must necessarily con- 
tain a local collective velocity term. Indeed the knowl- 
edge of the average spatial extension of the system at a 
given time, naturally produces a flow constraint at any 
successive time. Conversely a collective flow measure- 
ment at a given time can be translated into an informa- 
tion on the system density at a former time. 
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